PhytoREF: a reference database of the plastidial 16S rRNA gene of photosynthetic eukaryotes with curated taxonomy

Photosynthetic eukaryotes have a critical role as the main producers in most ecosystems of the biosphere. The ongoing environmental metabarcoding revolution opens the perspective for holistic ecosystems biological studies of these organisms, in particular the unicellular microalgae that often lack distinctive morphological characters and have complex life cycles. To interpret environmental sequences, metabarcoding necessarily relies on taxonomically curated databases containing reference sequences of the targeted gene (or barcode) from identified organisms. To date, no such reference framework exists for photosynthetic eukaryotes. In this study, we built the PhytoREF database that contains 6490 plastidial 16S rDNA reference sequences that originate from a large diversity of eukaryotes representing all known major photosynthetic lineages. We compiled 3333 amplicon sequences available from public databases and 879 sequences extracted from plastidial genomes, and generated 411 novel sequences from cultured marine microalgal strains belonging to different eukaryotic lineages. A total of 1867 environmental Sanger 16S rDNA sequences were also included in the database. Stringent quality filtering and a phylogeny-based taxonomic classification were applied for each 16S rDNA sequence. The database mainly focuses on marine microalgae, but sequences from land plants (representing half of the PhytoREF sequences) and freshwater taxa were also included to broaden the applicability of PhytoREF to different aquatic and terrestrial habitats. PhytoREF, accessible via a web interface (http://phytoref.fr), is a new resource in molecular ecology to foster the discovery, assessment and monitoring of the diversity of photosynthetic eukaryotes using high-throughput sequencing.

PhytoREF: a reference database of the plastidial 16S rRNA gene of photosynthetic eukaryotes with curated taxonomy

Photosynthetic eukaryotes have a critical role as the main producers in most ecosystems of the biosphere. The ongoing environmental metabarcoding revolution opens the perspective for holistic ecosystems biological studies of these organisms, in particular the unicellular microalgae that often lack distinctive morphological characters and have complex life cycles. To interpret environmental sequences, metabarcoding necessarily relies on taxonomically curated databases containing reference sequences of the targeted gene (or barcode) from identified organisms. To date, no such reference framework exists for photosynthetic eukaryotes. In this study, we built the PhytoREF database that contains 6490 plastidial 16S rDNA reference sequences that originate from a large diversity of eukaryotes representing all known major photosynthetic lineages. We compiled 3333 amplicon sequences available from public databases and 879 sequences extracted from plastidial genomes, and generated 411 novel sequences from cultured marine microalgal strains belonging to different eukaryotic lineages. A total of 1867 environmental Sanger 16S rDNA sequences were also included in the database. Stringent quality filtering and a phylogeny-based taxonomic classification were applied for each 16S rDNA sequence. The database mainly focuses on marine microalgae, but sequences from land plants (representing half of the PhytoREF sequences) and freshwater taxa were also included to broaden the applicability of PhytoREF to different aquatic and terrestrial habitats. PhytoREF, accessible via a web interface (http://phytoref.fr), is a new resource in molecular ecology to foster the discovery, assessment and monitoring of the diversity of photosynthetic eukaryotes using high-throughput sequencing.

Defining seasonal marine microbial community dynamics

Here we describe, the longest microbial time-series analyzed to date using high-resolution 16S rRNA tag pyrosequencing of samples taken monthly over 6 years at a temperate marine coastal site off Plymouth, UK. Data treatment effected the estimation of community richness over a 6-year period, whereby 8794 operational taxonomic units (OTUs) were identified using single-linkage preclustering and 21 130 OTUs were identified by denoising the data. The Alphaproteobacteria were the most abundant Class, and the most frequently recorded OTUs were members of the Rickettsiales (SAR 11) and Rhodobacteriales. This near-surface ocean bacterial community showed strong repeatable seasonal patterns, which were defined by winter peaks in diversity across all years. Environmental variables explained far more variation in seasonally predictable bacteria than did data on protists or metazoan biomass. Change in day length alone explains >65% of the variance in community diversity. The results suggested that seasonal changes in environmental variables are more important than trophic interactions. Interestingly, microbial association network analysis showed that correlations in abundance were stronger within bacterial taxa rather than between bacteria and eukaryotes, or between bacteria and environmental variables.

Ocean acidification shows negligible impacts on high-latitude bacterial community structure in coastal pelagic mesocosms

The impact of ocean acidification and carbonation on microbial community structure was assessed during a large-scale in situ costal pelagic mesocosm study, included as part of the EPOCA 2010 Arctic campaign. The mesocosm experiment included ambient conditions (fjord) and nine mesocosms with pCO(2) levels ranging from similar to 145 to similar to 1420 mu atm. Samples for the present study were collected at ten time points (t-1, t1, t5, t7, t12, t14, t18, t22, t26 to t28) in seven treatments (ambient fjord (similar to 145), 2x similar to 185, similar to 270, similar to 685, similar to 820, similar to 1050 mu atm) and were analysed for "small" and "large" size fraction microbial community composition using 16S rRNA (ribosomal ribonucleic acid) amplicon sequencing. This high-throughput sequencing analysis produced similar to 20 000 000 16S rRNA V4 reads, which comprised 7000OTUs. The main variables structuring these communities were sample origins (fjord or mesocosms) and the community size fraction (small or large size fraction). The community was significantly different between the unenclosed fjord water and enclosed mesocosms (both control and elevated CO2 treatments) after nutrients were added to the mesocosms, suggesting that the addition of nutrients is the primary driver of the change in mesocosm community structure. The relative importance of each structuring variable depended greatly on the time at which the community was sampled in relation to the phytoplankton bloom. The sampling strategy of separating the small and large size fraction was the second most important factor for community structure. When the small and large size fraction bacteria were analysed separately at different time points, the only taxon pCO(2) was found to significantly affect were the Gammaproteobacteria after nutrient addition. Finally, pCO(2) treatment was found to be significantly correlated (non-linear) with 15 rare taxa, most of which increased in abundance with higher CO2.

Minor impact of ocean acidification to the composition of the active microbial community in an Arctic sediment

Effects of ocean acidification on the composition of the active bacterial and archaeal community within Arctic surface sediment was analysed in detail using 16S rRNA 454 pyrosequencing. Intact sediment cores were collected and exposed to one of five different pCO(2) concentrations [380 (present day), 540, 750, 1120 and 3000 atm] and RNA extracted after a period of 14 days exposure. Measurements of diversity and multivariate similarity indicated very little difference between pCO(2) treatments. Only when the highest and lowest pCO(2) treatments were compared were significant differences evident, namely increases in the abundance of operational taxonomic units most closely related to the Halobacteria and differences to the presence/absence structure of the Planctomycetes. The relative abundance of members of the classes Planctomycetacia and Nitrospira increased with increasing pCO(2) concentration, indicating that these groups may be able to take advantage of changing pH or pCO(2) conditions. The modest response of the active microbial communities associated with these sediments may be due to the low and fluctuating pore-water pH already experienced by sediment microbes, a result of the pH buffering capacity of marine sediments, or due to currently unknown factors. Further research is required to fully understand the impact of elevated CO2 on sediment physicochemical parameters, biogeochemistry and microbial community dynamics.

Diverse profiles of N‐acyl‐homoserine lactone molecules found in cnidarians

Many marine habitats, such as the surface and tissues of marine invertebrates, including corals, harbour diverse populations of microorganisms, which are thought to play a role in the health of their hosts and influence mutualistic and competitive interactions. Investigating the presence and stability of quorum sensing (QS) in these ecosystems may shed light on the roles and control of these bacterial communities. Samples of 13 cnidarian species were screened for the presence and diversity of N-acyl-homoserine lactones (AHLs; a prevalent type of QS molecule) using thin-layer chromatography and an Agrobacterium tumefaciens NTL4 biosensor. Ten of 13 were found to harbour species-specific, conserved AHL profiles. AHLs were confirmed in Anemonia viridis using liquid chromatography tandem mass spectrometry. To assess temporal role and stability, AHLs were investigated in A. viridis from intertidal pools over 16 h. Patterns of AHLs showed conserved profiles except for two mid-chain length AHLs, which increased significantly over the day, peaking at 20:00, but had no correlation with pool chemistry. Denaturing gel electrophoresis of RT-PCR-amplified bacterial 16S rRNA showed the presence of an active bacterial community that changed in composition alongside AHL profiles and contained a number of bands that affiliate with known AHL-producing bacteria. Investigations into the quorum sensing-controlled, species-specific roles of these bacterial communities and how these regulatory circuits are influenced by the coral host and members of the bacterial community are imperative to expand our knowledge of these interactions with respect to the maintenance of coral health.

Response of an Arctic Sediment Nitrogen Cycling Community to Increased CO2

Ocean acidification influences sediment/water nitrogen fluxes, possibly by impacting on the microbial process of ammonia oxidation. To investigate this further, undisturbed sediment cores collected from Ny Alesund harbour (Svalbard) were incubated with seawater adjusted to CO2 concentrations of 380, 540, 760, 1,120 and 3,000 μatm. DNA and RNA were extracted from the sediment surface after 14 days' exposure and the abundance of bacterial and archaeal ammonia oxidising (amoA) genes and transcripts quantified using quantitative polymerase chain reaction. While there was no change to the abundance of bacterial amoA genes, an increase to 760 μatm pCO2 reduced the abundance of bacterial amoA transcripts by 65 %, and this was accompanied by a shift in the composition of the active community. In contrast, archaeal amoA gene and transcript abundance both doubled at 3,000 μatm, with an increase in species richness also apparent. This suggests that ammonia oxidising bacteria and archaea in marine sediments have different pH optima, and the impact of elevated CO2 on N cycling may be dependent on the relative abundances of these two major microbial groups. Further evidence of a shift in the balance of key N cycling groups was also evident: the abundance of nirS-type denitrifier transcripts decreased alongside bacterial amoA transcripts, indicating that NO3 − produced by bacterial nitrification fuelled denitrification. An increase in the abundance of Planctomycete-specific 16S rRNA, the vastmajority of which grouped with known anammox bacteria, was also apparent at 3,000 μatm pCO2. This could indicate a possible shift from coupled nitrification–denitrification to anammox activity at elevated CO2.

Single-cell enabled comparative genomics of a deep ocean SAR11 bathytype

Bacterioplankton of the SAR11 clade are the most abundant microorganisms in marine systems, usually representing 25% or more of the total bacterial cells in seawater worldwide. SAR11 is divided into subclades with distinct spatiotemporal distributions (ecotypes), some of which appear to be specific to deep water. Here we examine the genomic basis for deep ocean distribution of one SAR11 bathytype (depth-specific ecotype), subclade Ic. Four single-cell Ic genomes, with estimated completeness of 55%-86%, were isolated from 770 m at station ALOHA and compared with eight SAR11 surface genomes and metagenomic datasets. Subclade Ic genomes dominated metagenomic fragment recruitment below the euphotic zone. They had similar COG distributions, high local synteny and shared a large number (69%) of orthologous clusters with SAR11 surface genomes, yet were distinct at the 16S rRNA gene and amino-acid level, and formed a separate, monophyletic group in phylogenetic trees. Subclade Ic genomes were enriched in genes associated with membrane/cell wall/envelope biosynthesis and showed evidence of unique phage defenses. The majority of subclade Ic-specfic genes were hypothetical, and some were highly abundant in deep ocean metagenomic data, potentially masking mechanisms for niche differentiation. However, the evidence suggests these organisms have a similar metabolism to their surface counterparts, and that subclade Ic adaptations to the deep ocean do not involve large variations in gene content, but rather more subtle differences previously observed deep ocean genomic data, like preferential amino-acid substitutions, larger coding regions among SAR11 clade orthologs, larger intergenic regions and larger estimated average genome size.

Rapid response of the active microbial community to CO2 exposure from a controlled sub-seabed CO2 leak in Ardmucknish Bay (Oban, Scotland)

The response of the benthic microbial community to a controlled sub-seabed CO2 leak was assessed using quantitative PCR measurements of benthic bacterial, archaeal and cyanobacteria/chloroplast 16S rRNA genes. Samples were taken from four zones (epicentre; 25 m distant, 75 m distant and 450 m distant) during 6 time points (7 days before CO2 exposure, after 14 and 36 days of CO2 release, and 6, 20 and 90 days after the CO2 release had ended). Changes to the active community of microphytobenthos and bacteria were also assessed before, during and after CO2 release. Increases in the abundance of microbial 16S rRNA were detected after 14 days of CO2 release and at a distance of 25 m from the epicentre. CO2 related changes to the relative abundance of both major and minor bacterial taxa were detected: most notably an increase in the relative abundance of the Planctomycetacia after 14 days of CO2 release. Also evident was a decrease in the abundance of microbial 16S rRNA genes at the leak epicentre during the initial recovery phase: this coincided with the highest measurements of DIC within the sediment, but may be related to the release of potentially toxic metals at this time point.

Metabolically active, non-nitrogen fixing,Trichodesmiumin UK coastal waters during winter

Trichodesmium, a colonial cyanobacterium typically associated with tropical waters, was observed between January and April 2014 in the western English Channel. Sequencing of the heterocyst differentiation (hetR) and 16S rRNA genes placed this community within the Clade IV Trichodesmium, an understudied clade previously found only in low numbers in warmer waters. Nitrogen fixation was not detected although measurable rates of nitrate uptake and carbon fixation were observed. Trichodesmium RuBisCO transcript abundance relative to gene abundance suggests the potential for viable and potentially active Trichodesmium carbon fixation. Observations of Trichodesmium when coupled with a numerical advection model indicate that Trichodesmium communities can remain viable for >3.5 months at temperatures lower than previously expected. The results suggest that Clade IV Trichodesmium occupies a different niche to other Trichodesmium species, and is a cold- or low-light-adapted variant.

Dynamic responses of the benthic bacterial community at the Western English Channel observatory site L4 are driven by deposition of fresh phytodetritus

The impact of the seasonal deposition of phytoplankton and phytodetritus on surface sediment bacterial abundance and community composition was investigated at the Western English Channel site L4. Sediment and water samples were collected from January to September in 2012, increasing in frequency during periods of high water column phytoplankton abundance. Compared to the past two decades, the spring bloom in 2012 was both unusually long in duration and contained higher than average biomass. Within spring months, the phytoplankton bloom was well mixed through the water column and showed accumulations near the sea bed, as evidenced by flow cytometry measurements of nanoeukaryotes, water column chlorophyll a and the appearance of pelagic phytoplankton at the sediment. Measurements of chlorophyll and chlorophyll degradation products indicated phytoplankton material was heavily degraded after it reached the sediment surface: the nature of the chlorophyll degradation products (predominantly pheophorbide, pyropheophorbide and hydroxychlorophyllone) was indicative of grazing activity. The abundance of bacterial 16S rRNA genes g−1 sediment (used as a proxy for bacterial biomass) increased markedly with the onset of the phytoplankton bloom, and correlated with measurements of chlorophyll at the surface sediment. Together, this suggests that bacteria may have responded to nutrients released via grazing activity. In depth sequencing of the 16S rRNA genes indicated that the composition of the bacterial community shifted rapidly through-out the prolonged spring bloom period. This was primarily due to an increase in the relative sequence abundance of Flavobacteria.

Seasonal variability in microbial methanol utilisation in coastal waters of the western English Channel

Methanol is ubiquitous in seawater and the most abundant oxygenated volatile organic compound (OVOC) in the atmosphere where it influences oxidising capacity and ozone formation. Marine methylotrophic bacteria utilise methanol in seawater both as an energy and/or growth substrate. This work represents the first fully resolved seasonal study of marine microbial methanol uptake dynamics. Rates of microbial methanol dissimilation in coastal surface waters of the UK varied between 0.7 – 11.2 nmol l-1 h-1 and reached a maximum in February. Rates of microbial methanol assimilation varied between 0.04 – 2.64 x 10-2 nmol l-1 h-1 and reached a maximum in August. Temporal variability in microbial methanol uptake rates shows that methanol assimilation and dissimilation display opposing seasonal cycles, although overall <1% of methanol was assimilated. Correlative approaches with 16S rRNA pyrosequencing data suggested that bacteria of the SAR11 clade and Rhodobacterales could be significantly influencing rates of methanol dissimilation and assimilation, respectively, at station L4 in the western English Channel

Seasonal variability in microbial methanol utilisation in coastal waters of the western English Channel

Methanol is ubiquitous in seawater and the most abundant oxygenated volatile organic compound (OVOC) in the atmosphere where it influences oxidising capacity and ozone formation. Marine methylotrophic bacteria utilise methanol in seawater both as an energy and/or growth substrate. This work represents the first fully resolved seasonal study of marine microbial methanol uptake dynamics. Rates of microbial methanol dissimilation in coastal surface waters of the UK varied between 0.7 – 11.2 nmol l-1 h-1 and reached a maximum in February. Rates of microbial methanol assimilation varied between 0.04 – 2.64 x 10-2 nmol l-1 h-1 and reached a maximum in August. Temporal variability in microbial methanol uptake rates shows that methanol assimilation and dissimilation display opposing seasonal cycles, although overall <1% of methanol was assimilated. Correlative approaches with 16S rRNA pyrosequencing data suggested that bacteria of the SAR11 clade and Rhodobacterales could be significantly influencing rates of methanol dissimilation and assimilation, respectively, at station L4 in the western English Channel